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Abstract 

The Lifshitz theory of the van der Waals force is extended for the case of an atom (molecule) 
interacting with a plane surface of an uniaxial crystal or with a long solid cylinder or cylindrical 
shell made of isotropic material or uniaxial crystal. For a microparticle near a semispace or flat 
plate made of an uniaxial crystal the exact expressions for the free energy of the van der Waals 
and Casimir-Polder interaction are presented. An approximate expression for the free energy 
of microparticle-cylinder interaction is obtained which becomes precise for microparticle-cylinder 
separations much smaller than cylinder radius. The obtained expressions are used to investigate 
the van der Waals interaction between hydrogen atoms (molecules) and graphite plates or multiwall 
carbon nanotubes. To accomplish this the behavior of graphite dielectric permittivities along the 
imaginary frequency axis is found using the optical data for the complex refractive index of graphite 
for the ordinary and extraordinary rays. It is shown that the position of hydrogen atoms inside 
multiwall carbon nanotubes is energetically preferable compared with outside. 
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I. INTRODUCTION 



The van der Waals interaction between microparticle and macrobody has long been inves- 
tigated. It is of much importance for understanding of a large body of physical and chemical 
phenomena connected with atom-surface interaction including adsorption and friction. In a 
pioneering work in Ref. the interaction potential between an atom at a separation a from 
a plane wall was found in the from V^la) = —C^/a^. This result is applicable at separations 
less than a few nanometers. More recently, a lot of different atoms, molecules and wall 
materials was studied. In particular, in Refs. the values of C3 were computed for the 

interaction of H, H2, He, Ne, Ar, Cr, Xe, and CII4 with the planar surfaces of insulators (sap- 
phire, LiF, CaF2, and boron nitride). At much greater separations the atom- wall interaction 
is described by the Casimir-Polder potential V^la) = —Ci/a^ 0] taking relativistic effects 
into account. The complete theory of the van der Waals atom- wall interaction at nonzero 
temperature is given by the Lifshitz formula in terms of the dynamic polarizability of an 
atom (molecule) and the frequency- dependent dielectric permittivity of wall material. The 
potentials V^{a) and V4(a), obtained previously, are the two limiting cases of this formula. 

During the last few years van der Waals forces have found important new applications in 



experiments on quantum reflection and diffraction o 



and in Bose-Einstein condensation 
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ultra-cold atoms on different surfaces 



In connection with this, the detailed 
examination of different corrections to the Casimir-Polder and van der Waals interactions, 
including the precise effect of atomic polarizability and nonideality of wall material was per- 
formed in Refs. 12 , ^ . Effectively this resulted in the investigation of accurate dependences 
of the coefficients C3 and C4 on separation and temperature. 

Although the Lifshitz theory presents considerable opportunity for extensive studies of 



the van der Waals force 



15l |. it is essentially restricted by macroscopic bodies with plane 



boundaries. The use of approximations, like the proximity force theorem |16j], permitted 
one to obtain rather precise results for a large sphere near a plane plate, a confi gura tion 
frequently used in recent experiments on measuring the Casimir force Q Q Q H, Q- 
In most cases the macrobodies with plane boundaries were supposed to be isotropic. 

In the present paper we generalize the Lifshitz formula for a microparticle situated near 
the surface of an uniaxial crystal. Both cases of crystal semispace with plane boundary and 
a plane plate of finite thickness are considered. As a next step, we derive the approximate 
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expression for the free energy of the van der Waals interaction between a microparticle and 
a sohd cyhnder or cyhndrical shell made of an uniaxial crystal. In the limiting case this 
expression is applicable to a microparticle near a cylinder made of an isotropic material 
with frequency dependent dielectric permittivity (a configuration which also has not been 
investigated previously). We apply the obtained results to investigate the van der Waals 
interaction between hydrogen atoms or molecules and graphite plates or multiwall carbon 
nanotubes. 

The study of the van der Waals interaction between hydrogen atoms and a graphitic sur- 
face has become urgent after the proposal of Ref. [3] to use the singlewall carbon nanotubes 
for hydrogen storage. Since, many papers were published on the use of both singlewall and 
multiwall nanotubes for hydrogen storage and containing both promising and disappointing 
results (see Ref. ||23| for review). The macroscopic theoretical approach leads to a conclusion 



2J| that the carbon nanostructures might absorb hydrogen from 4 to 14 percent of their 



weight. However, the microscopic mechanisms responsible for this absorption are still un- 
known. The van der Waals forces acting between hydrogen atoms or molecules and carbon 
nanostructures, which might play an important role in absorption phenomena, are practi- 



cally unexplored. Some pre 



be found in Refs. 23 and 
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iminary results for graphite sheets and singlewall nanotubes can 
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27|, respectively. The van der Waals interaction of fulerene 



molecules and adsorption of these molecules on graphite were considered in Ref. j28|. 

To apply the Lifshitz-type formulas for the van der Waals free energy, obtained in the 
paper, to the case of hydrogen atoms and molecules near graphite surface, we calculate 
the dielectric permittivities of graphite and dynamic polarizabilities of hydrogen atom and 
molecule along the imaginary frequency axis. To do this, we discuss different sets of tabulated 
optical data for the complex refractive index of graphite and use the most reliable ones to 
perform the Kramers-Kronig analysis. The van der Waals interactions between hydrogen 
atom and molecule and graphite semispace or plate of finite thickness are calculated. The 
free energies of hydrogen atom inside and outside of a multiwall carbon nanotube are found 
as functions of an atom-nanotube separation distance and internal and external nanotube 
radia. The location of a hydrogen atom inside a multiwall nanotube is demonstrated to be 
preferable from an energetic point of view. 

The paper is organized as follows. In Sec. II we present the Lifshitz formula for the van 
der Waals (and Casimir-Polder) interaction between microparticle and plane surface of an 
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uniaxial crystal. Sec. Ill contains derivation of general expression for the van der Waals free 
energy of a microparticle external to a solid cylinder or cylindrical shell made of an uniaxial 
crystal. In Sec. IV the dielectric permittivities of graphite and the atomic and molecular 
dynamic polarizabilities of hydrogen along the imaginary frequency axis are obtained. In 
Sec. V calculation results are presented for the van der Waals interaction between hydrogen 
atom or molecule and graphite semispace or a plane plate of finite thickness. In Sec. VI the 
same is done for hydrogen atom or molecule external to a multiwall carbon nanotube. Com- 
parison between the free energies of hydrogen atom inside and outside multiwall nanotube 
is done in Sec. VII. Sec. VIII contains our discussion and conclusions. 



II. LIFSHITZ FORMULA FOR THE VAN DER WAALS INTERACTION 
BETWEEN MICROPARTICLE AND PLANE SURFACE OF AN UNIAXIAL 
CRYSTAL 

First we consider a neutral microparticle (atom or molecule) with a dynamic polariz- 
ability a{uj) at separation a from a plane surface of the isotropic semispace with dielectric 
permittivity e{uj) at temperature T in thermal equilibrium. In this case the free energy of 
microparticle-semispace van der Waals interaction is given by the familiar Lifshitz formula 
B (see also 11, 12, 3, HQ) 



-2aqi 



F^ia^T) = -ksT } ^ a{iC,i) / k±dk±qie~ 
1=0 ^0 

X |2r^(e^ k^) + -|_ [rli^^^ k^) _ rl{iu k^)] } • (1) 

Here = 2iTkBTl/h are the Matsubara frequencies, ks is the Boltzmann constant, / = 
0, 1, 2, . . . , and k± is the magnitude of a wave vector component in the plane surface of a 
semispace. The coefficients of reflection for two independent polarizations of electromagnetic 
field are given by 

£iqi - h 



rUuk, 



(2) 



4 



where 



£i 



e + 



kl + €i 



(3) 



and prime near the summation sign in Eq. (P) means that the term for I = has to be 
multiphed by 1/2. 

Eq. can be readily generahzed for the case when the microparticle is located not 
near a semispace, but near a flat plate of some finite thickness d with the same dielectric 
permittivity e{uj). In this case the free energy of the van der Waals interaction F^{a, T) again 
is given by Eq. (0) where, however, the reflection coefficients from a semispace ^{^i^k^^ 
should be replaced by the reflection coefficients from a plate of finite thickness '"[i k±). 
The explicit expressions for them are obtained from the free energy of the van der Waals 
interaction between the layered media (see, e.g.. 



efqf + kf + 2qikieicoth.{kid) ' 



kf - qf 



(4) 



kf + qf + 2qikiCoth.{kid) ' 
In the limit c/ — > oo Eq. (0} transforms into Eq. (0)- 

Let us now consider a semispace or a plate of finite thickness made of an uniaxial crys- 
tal (graphite for instance) which is characterized by two dissimilar dielectric permittivities 
Exijjj) = ey{u) and ez{uj). Let a microparticle be located near the uniaxial crystal semispace 
restricted by the plane {x,y), and the crystal optical axis z being perpendicular to it. Then 
the free energy of the van der Waals interaction is again given by Eq. (P) where the coeffi- 
cients of reflection from the surface of isotropic semispace ry j_(^k±) should be replaced by 
their generalization for the case of uniaxial crystal (graphite 

y/^xl^zlQl — kzi 



rl,gi^l,k±) : 

rl-Muk^) 



Here the following notations are introduced 



y/^xl^zlQl + kzi 

kxi - qi 
kxi + qi ' 



(5) 



kxi 

£xl 



A;2 +£ 



xl ' 



£x{i^i), 



kzi = 
£z{i^i)- 



k^, +e 



zl- 



(6) 



For isotropic crystal Ex = Sz = e and Eq. (0) coincides with Eq. (j21). 

If a microparticle is located near a flat plate of finite thickness made of uniaxial crystal 
(z-axis is perpendicular to the plate), the free energy F^{a,T) is given again by Eq. (P), 
where the coefficients of reflection from an isotropic plate r^^ kj_) are replaced by the 
reflection coefficients from a plate made of uniaxial crystal: 



SxiBzigf + + 2y/6xieziqikzicoth{kzid) ' 

kli + qf + 2qikxiCothikxidy ^' ^ 

For the anisotropic plate of infinite thickness {d —>■ oo) Eq. ((7j) transforms into Eq. (jS)). On 
the other hand, in the limit of the plate made of isotropic substance Eq. ((7j) coincides with 
Eq. dH). 

Eq. ((H) with reflection coefficients ©, ((Zj) is used in Sec.V for computations of the van 
der Waals interaction between the hydrogen atoms or molecules and the plane surface of a 
semispace or a plate made of graphite. 

III. FREE ENERGY OF THE VAN DER WAALS INTERACTION FOR 
A MICROPARTICLE EXTERNAL TO A SOLID OR HOLLOW CYLINDER 

In this section we derive the Lifshitz-type formula for the van der Waals free energy of 
a microparticle located at a separation a from the external surface of a solid cylinder or 
cylindrical shell made of an uniaxial crystal. It is assumed that the crystal optical axis z 
is perpendicular to the cylinder surface of crystalline layers. The outer radius of a cylinder 
is R and the thickness of a crystal cylindrical shell is d < R. In the case d = R the 
cylinder is solid. If d < R, there is an empty cylindrical cavity inside of a cylinder. As 
in the previous section, the crystalline material of the cylindrical shell is described by the 
dielectric permittivities e^iuj) and e^iuj). The derivation presented below is based on the 
same approach which was previously used in literature p, 12, 29, 30, [sil to derive the 
Lifshitz formula for microparticle-semispace (plate) interaction from the Lifshitz formula for 
a configuration of two parallel semispaces (plates). 

Let us consider an infinite space filled with an isotropic substance having a dielectric 
permittivity e{uj), containing an empty cylindrical cavity of radius R + a. We introduce 
our solid cylinder or cylindrical shell of external radius R made of an uniaxial crystal inside 



this cavity so that the cyhnder axis coincides with the axis of the cavity (see Fig. 1). Then 
there is a gap of thickness a between our cyhnder and the boundary of the cyhndrical cavity 
of radius R + a restricting the infinite space with the dielectric permittivity s{uj). Each 
element of our cylinder experiences an attractive van der Waals interaction on the source 
side of the boundary of the cylindrical cavity restricting the infinite space. With the help 
of the proximity force theorem the free energy of this interaction between two cylinders can 
be approximately represented in the form (see Ref. js^ for the case of ideal metals) 



F'gXa, T) = 2nLy/R{R + a)F'/{a, T). (8) 

Here F^'^la, T) is the free energy per unit area in the configuration either of two semispaces 
separated by a gap of width a (in this case i = s, our cylinder is solid, one semispace is filled 
with an uniaxial crystal and the other is filled with a material of dielectric permittivity e{u!)) 
or of a flat plate of thickness d and a semispace separated by the same gap (in this case 
i = p, and we are dealing with cylindrical shell having a longitudinal hole of radius R — d; the 
plate is made of an uniaxial crystal and semispace of material with a dielectric permittivity 
e{uj)). In Eq. (jHJ L is the length of our solid or hollow cylinder which is supposed to be 
much larger than its radius R. 

As shown in Ref. 1351 (see also Ref. 



the accuracy of Eq. (jSI) is rather high. For exam- 
ple, within the separation region < a < R/2 the results calculated by Eq. (jH)) coincide with 
the exact ones up to 1% in the case of cylinders made of perfect metal (for other materials 
the accuracy may be different for only a fraction of percent). This is quite satisfactory for 
application to multiwall nanotubes with R of about a few ten nanometers considered below. 
The explicit expressions for the free energy F^£{a, T) are well known 

F^^{a,T) = ^Y. / ^A{in[i-^;:^te,fc±)i(e/,fc±)e-''^^' 

+ In [1 - rl:f^(ez, k^YAiu k^)e'^'''^^] ] . (9) 

Here the reflection coefficients r^^ ^.^ from the semispace of uniaxial crystal are given by 
Eq. (0), coefficients ±.g, describing reflection from a flat plate of uniaxial crystal, are given 
by Eq. ((7j), and coefficients r^^ ^ describing reflection from isotropic semispace are presented 
in Eq. Q. Notice that when index i in the left-hand side of Eq. Q is equal to s or p one 
should choose s or p in the right-hand side, respectively. 
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To continue with our derivation, we now suppose that the isotropic substance with the 
dielectric permittivity e{u!) is rarefied with the number of atoms or molecules per unit 
volume. Expanding the quantity -F^'^(ci, T) from the left-hand side of Eq. (jSj) as a power 
series in and using the additivity of the first-order term, one can write 



POO 

F'/{a, T) = N F^{z, T)2n{R + z)Ldz + 0{N^' 

J a 



(10) 



where F§;{z, T) is the free energy of the van der Waals interaction of a single atom belonging 
to an isotropic substance with a solid cylinder or cylindrical shell made of an uniaxial crystal 
(note that separation z is measured from the external surface of the cylinder in the direction 
perpendicular to it). 

By differentiation of both sides of Eq. (fTUI) with respect to a, we obtain 

dF^/{a,T) 



da 



2tt{R + a)LNF^{a, T) + 0{N^ 



(11) 



The same derivative can be found when differentiating both sides of Eq. 

9F|^(a,T) 



da 



2ti L^R{R + 
1 



(12) 



2{R + a] 



where 



F'''{a,T) 



F'i\a,T) + F^'%a,T) 



dFj^%a,T) 



da 



(13) 



is the van der Waals force per unit area acting between the semispace made of an uniaxial 
crystal {i = s) or a flat plate made of the same material and a semispace with a dielectric 
permittivity e. The expression for this force is easily obtained from Eqs. Q and (fT^ : 



F^'^(a,T) = -^V / k^dk^qi 



(14) 



X 



1=0 



+ 



rZ{^i,k±)rl{^i,k^) 



The dielectric permittivity of a rarefled substance can be expanded in Taylor series in 
powers of js^] 

em = 1 + Ana{i^i)N + 0{N^), (15) 
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where a{uj) is the dynamic polarizabihty of an atom (molecule) of this substance. Substi- 
tuting Eq. (fT3j) in Eqs. ^ and ^ we obtain 



qfc^ 



+ 0{N'^ 



(16) 



Using Eq. the free energy F^'^ and the force F*''' from Eqs. and (jl4p can be 

represented in the form 



F^^(a,T) 



k\dk\ 



qfc'^ 



'2aqi 



+ 0{N'^ 



F^'^(a,T) = -fcsTiVE a(2eO / k^dk^qi 

1=0 -^0 



(17) 



2- J^) + 



+ 0(A^2 



Substituting Eq. (HH) in Eq. jHl), one finds 
9F|^(a,T) _ 



5a 



2^TLNkBT^/R{R + a) ^ «(«^i) / 



A{R + a) 



X <; 2r^;^(e, fc^) + ri:f^(e. fc^) - r;:^(e^ 



qfc'^ 



(18) 



As a final stage of the derivation, we substitute the result (jl8p into the left-hand side of 
Eq. take the limit N ^ and arrive at desired expression for the free energy of van 
der Waals interaction between a microparticle and a cylinder made of uniaxial crystal 



Flj{a,T) = -ksT 



R 



R + a 



oo ^ „c 
1=0 



A{R + a) 



X <! 2r;;:^(6, fc^) + ^ rl:f^(6, A;x) - rZi^i, k 



qfc 



(19) 



In the case of a solid cylinder, the reflection coefficients f\\±.g, given by Eq. (jSl), should be 



chosen in the right-hand side of Eq. ()19|). For a cylindrical shell, coefficients ^^ii j^.^ from 



Eq. (|7j) should be used. Notice that in the limit R ^ oo Eq. (fT^ coincides with a known 
result for the free energy of microparticle near a plane surface of a semispace. The 
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above derivation is preserved also in the limiting case of a solid or hollow cylinder made of 
isotropic material with Ex = Ey = Sz = s. To obtain the result for isotropic cylinder, one 
should substitute in Eq. (fT^ the reflection coefficients 0, dH) instead of (0), (0). 

Eq. (jl9p is the approximate one. It is, however, practically exact at a ^ i? and is of 
high precision (the error is of about 1%) at all separations a < R/2. That is why this 
equation is reliable for calculations of the van der Waals interaction between a cylinder and 
microparticles located in its close proximity. 



IV. DIELECTRIC PERMITTIVITIES OF GRAPHITE AND DYNAMIC 
POLARIZABILITIES OF HYDROGEN ATOM AND MOLECULE ALONG THE 
IMAGINARY FREQUENCY AXIS 

Below we used the Lifshitz-type formulas obtained above to calculate the van der Waals in- 
teraction between hydrogen atoms or molecules and graphite semispace or flat plate [Eqs. (^J), 
®, ©] or graphite cyhnder [Eqs. ©, (j?!), (fT^ ]. The graphite cylinder models a multiwall 
carbon nanotube (see Sec. VI). To attain these ends, one needs the values of dynamic polar- 
izabilities of hydrogen atom and molecule and also both dielectric permittivities of graphite 
at all Matsubara frequencies which give non- negligible contribution to the result. 

The precise expression for the atomic dynamic polarizability of hydrogen is given by the 



10-oscillator formula 



38| written in atomic units 
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«(^^') = E-7^' (20) 

j=l "^aj + 4z 

where Qj are the oscillator strengths and uJaj are the eigenfrequencies. For the hydrogen atom 
the values of these quantities are listed in Table I (note that 1 a.u. of energy = 4.3597 x 
10^^^ J = 27.11 eV). Note also that before the substitution in Eqs. (0) or ()19p the atomic 
dynamic polarizability from Eq. (j^UI) should be expressed in cubic meters including the 
transformation factor for 1 a.u. of polarizability = 1.482 x 10~^^ m^. 

In addition to the precise representation ()2()|1 . the atomic dynamic polarizability of hy- 
drogen atom can be expressed in terms of a more simple single oscillator model 

= -TTT-^ (21) 
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where ga = aa(0)ci;^ is expressed through the static atomic polarizabihty aa{0) = 4.50 a. u. 
and the characteristic energy uJa = 11.65 eV 39^ . 

Below we will check that after the substitution to the Lifshitz-type formulas both expres- 
sions (j2(J|) and (j2H) lead to equal results in the limits of required accuracy. This permits to 
use a more simple Eq. (j2ip in computations. 

It is well known that for hydrogen molecule the single oscillator model for the dynamic 
polarizabihty is more exact than for the atom. For this reason it is acceptable to present 
the molecular dynamic polarizabihty of hydrogen in the form 

where gm = o:m{0)u!'^. Here the static polarizabihty and the characteristic energy of hydro- 
gen molecule are equal to 0^(0) = 5.439 a.u. and lu^ = 14.09 eV, respectively 39|. 

Now let us consider the problem of dielectric permittivities of graphite Ex and Ez along 
the imaginary frequency axis. Both these quantities can be computed with the help of 
Kramers-Kronig relation 

^ Jo + ? 

The imaginary parts of the respective dielectric permittivities along the real axis, in turn, 
are equal to 2Iienx,z{^) x i-e., are expressed through the real and imaginary parts 

of the c omp lex refractive index of graphite for ordinary and extraordinary rays, respectively. 

Ref. |40| contains the measurement data for Kenx^zi^^) and Imux^zi^) of graphite obtained 
by different authors in the frequency region from Qi = 0.02 eV to ^2 = 40 eV (1 eV = 1.519 x 
lO^^rad/s). The use of these data to calculate Ex,z{iO by Eq- is, however, complicated 
by the two problems. First, the interval [fii,i72] is too narrow to calculate E^^zi'^O 
Matsubara frequencies contributing to the van der Waals force (by comparison, for Au the 
complex refractive index is measured from 0.125 eV to 10000 eV). Second, although for 
data by different authors are in agreement, in the case of Uz there are contradictory data in 
literature at u; < 15.5 eV. 

The first problem can be solved by the use of extrapolation. According to Ref. Q], at 
high frequencies u > Q2 the imaginary parts of graphite dielectric permittivities can be 
presented analytically in the form 

I^^S(^) = (24) 



a;3 
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Here the values of constants Ax = 9.60 x 10^ eV^ and = 3.49 x 10^ eV^ are determined 
from the condition of a smooth joining with the tabulated data at u = Q2 13 • 

At low frequencies u < Qi one may approximate Ime^ with the help of the Drude model 



Ime«(, 



(25) 
0.04 eV are 



uj^uj"^ + 7^) 

where the plasma frequency uOp = 1.226 eV and the relaxation parameter 7 
determined from the demand of smooth joining with tabulated data at u = Qi. 

The extrapolation of tabulated data for Ime^ to the region of low frequencies is connected 
with the second problem discussed above, i.e., with the contradictory measurements by 
different authors. Thus, the measurement data for 11^(00) in Ref. {4^ differ considerably 



from the same data in Ref. 
Refs 



41 



42| in the frequency region u < 15.5 eV. According to both 



42| |. the imaginary part of ^^(ci;) can be extrapolated to low frequencies uj < Qi 
by a constant: 



Imel'^i 



'zO 



const. 



The values of this constant, however, are found to be different: e'!.Q 



(26) 

3 according to Ref. [4 



and e'!.Q = according to Ref. |41|. 

As a result, the calculation of graphite dielectric permittivities along the imaginary fre- 
quency axis by Eq. (pHj) is performed as follows: 



rc .In UJ^ + 4^ 



^2 

Hi 



^0 

cjlmei^l 



, culmeij 
n2 ^ + i 



where Ime^*^^ is found from the tables and Ime^x^'P are given by Eqs. 

Eqs. (I21-(j2ni) in Eq. ^ one finds 

2 4'Arctan^ - 7Arctan^ 
ex{ii) = 1 + TT^TT^ -UJ^^ 



(27) 

Substituting 



TT 



+ 



duj 



uj\m.ex\i^) 



7 



+ 



Ax 



2 1 

+ 



^.(^0 = l + ^lnfl + ^ 

V ^ 



no 



2 , ^2 
Arctan 1 



(28) 



H 



21/2 ^2 

—— + - -Arctan— 



1 
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The calculational results from Eq. 



obtained by the use of the tabulated optical data 
of Refs. |41|,|42], are shown in Figs. 2a,b in the frequency range from = 2.47x 10^^ rad/s 
to ^2000 at T = 300K. These results allow the precise calculation of the van der Waals 
interaction by Eqs. (Q), (jl9p in the separation region a > 3 nm (note that with the increase 
of separation the number of Matsubara frequencies, giving a non-negligible contribution to 
the result, decreases). As to the contribution of zero Matsubara frequency = 0, there is the 
analytical result '"i|.g(0, k±) = 1 which follows from e^l^O ~^ ^ when ^ ^ in accordance 
with Eq. (|^. Note that at zero frequency the other reflection coefficient r^''^(0, does 
not contribute to the result due to the multiple the right-hand sides of Eqs. ((H) and 

(Hi). 

The de pen dence of e^iiO ^ Fig- 2a is typical for good conductors (compare with 
33| for Al and Au). In Fig. 2b the solid line is obtained with the results of 
see also Ref. 




4ll (see also Ref. 



Ad ) with e'!,Q = 3. The dashed line in Fig. 2b is obtained by the data 
Q|) using e'^Q = 0. It is seen that the dashed line differs markedly 
from the solid line in the frequency region ^ < 10^'' rad/s. The respective differences in the 
free energy are discussed in the next section. It is reasonably safe, however, to prefer the 
solid line in Fig. 2b as giving the correct behavior of along the imaginary frequency axis. 
In fact the difference between the two lines is due to the absence of absorption bands near 



the frequencies of 5eV and 11 eV in the tabulated data of Ref. 



4l| related to (note that 



in the data for Ex there are absorption bands at these frequencies in both Refs. |4ll . 
This casts doubts on the measurement data of Ref. for Ez because from the theory of 
graphite band structure j^l it follows that the respective absorption bands must be present 
simultaneously in both sets of data for e^ and Ez- 



V. CALCULATION OF THE VAN DER WAALS INTERACTION BETWEEN 
HYDROGEN ATOM OR MOLECULE AND PLANE SURFACE OF GRAPHITE 

We consider the hydrogen atom or molecule at a separation a from the hexagonal plane 
surface (x, y) of a graphite semispace of a fiat graphite plate of thickness d. Note that the 
separation distance between the two plane hexagonal layers in graphite is approximately 
0.336 nm. All calculations are performed at separations a > 3nm where one can neglect the 
atomic structure of graphite and describe it in terms of dielectric permittivities Ex{uj), Ez{uj) 
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as is done in the Lifshitz theory. Bearing in mind apphcations at short separations, it is 
instructive to present Eq. in the form of nonrelativistic van der Waals interaction (see 
Introduction) 

F.P(a,T) = -S!(^, (29) 

Qj 

where the van der Waals coefficient C|'^ [for the case of an atom near a semispace (s) or a 
plate (p), respectively] is now a function of both separation and temperature. For the sake 
of convenience in numerical computations, we introduce the nondimensional variables 

y = 2aqu = ^ = ^ (30) 

C UJc 

and express the van der Waals coefficient in terms of these variables 

^'^(a, T) = ^l 2a(0) + ^ a{iQu:,) (31) 
I 1=1 



CO 

y 



X / dye 



22/V;:^(0,2/) + cf ri:f^(0,2/)-r;;:^(0,2/) 



Note that for separations up to a few hundred nanometers Eq. (PT|) practically does not 
depend on temperature. 

In terms of the new variables (j3(jp the coefficients of reflection from a graphite semispace 
(0) are rearranged as 

^/£xi£ziy - fzivXi) 



rL{Ci,y) 



(32) 

fx{y,Ci) - y 
fx{yXi) + y' 

where 

/'(2/,0) = y' + Cf(^./-i), 

/'(2/,0)=y' + Cf(^./-l). (33) 
In analogy, the reflection coefficients (|7j) from a flat plate of thickness d take the form 

p /> N ^xi^ziy'^ — fziyXi) 

"""^^^ " SxiSziy^ + Pz{y. Ci) + '^V^iyf^^y^ C/)coth [/,(y, Ci)d/{2a)] ' 



(34) 

Px{y.Q)-y' 



y' + P{y. 0) + 2i//.(i/, O)coth [/^y, 0)rf/(2a)] ' 
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Now we substitute the reflection coefficients from a semispace (j52jl . tlie precise atomic 
dynamic polarizability ()20|) and data of Fig. 2a for Sx and Fig. 2b (solid line) for Sz into 
Eq. (j3ip . The calculational results for the coefficient of van der Waals interaction between 
a hydrogen atom and graphite semispace are presented in Fig. 3a by the solid line. For 
comparison the dashed line in Fig. 3a shows the results obtained with the use of alternative 
data for Sz (dashed line in Fig. 2b). As is seen from Fig. 3a, at the shortest separation 
a = 3 nm the use of the alternative data for Ez leads to a 15% error in the value of the van 
der Waals coefficient which decreases with an increase of separation. 

The computation of C| was repeated using the single oscillator model ^2\\ for the atomic 
dynamic polarizability instead of the 10-oscillator model (|2(Jj). The results were found to be 
practically in coincidence with those in Fig. 3a (the maximum deviations are less than 0.2% 
in the separation region from 3nm to 150 nm). Thus, the single oscillator model is a suf- 
ficient approximation for the atomic (and, consequently, molecular) dynamic polarizability 
of hydrogen in computations of the short-range van der Waals interaction with a graphite 
surface. 

In the same way as above, we calculate the van der Waals coefficient C| for the interaction 
of a hydrogen molecule with graphite semispace. The only difference is the use of the 
molecule dynamic polarizability (j^^ instead of atomic one. The results are shown in Fig. 3b 
by the solid line (the dashed line is calculated by the less accurate alternative data of Ref. 41 1 
for the dielectric permittivity £z)- The comparison of Figs. 3a and 3b leads to the conclusion 
that the magnitudes of the van der Waals coefficient for the hydrogen molecule are larger 
than for the atom. 

Now let the hydrogen atom be located at a separation a from the flat graphite plate of 
thickness d. Of interest is the dependence on d of the van der Waals free energy of atom- 
plate interaction. The calculations of the free energy were performed by Eqs. (j^^ and ()31|) 
with reflection coefficients (j32|) (for a semispace) and p4|l (for a plate of thickness d). The 
values of dielectric permittivities along the imaginary frequency axis were taken from Fig. 2 
(solid lines) and the atomic dynamic polarizability from Eq. ()2H) . In Fig. 4 the ratios of 
the free energies are plotted for the case of a plate and a semispace as a function of plate 
thickness for hydrogen atom located at different separations from the graphite surface (line 
1 for a = 3nm, line 2 for a = 10 nm, line 3 for a = 20 nm, and line 4 for a = 50 nm). 
As is seen from Fig. 4, at a separation a = 3nm the flnite thickness of the plate has a 
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pronounced effect on the free energy (more than 1% change) only for thcknesses d < 8nm. 
At separations a = 10 nm, 20 nm and 50 nm the finite thickness of the plate leads to a smaller 
magnitude of the van der Waals free energy, as compared with a semispace, for more than 
1% if the thickness of a plate is less than 19 nm, 32 nm and 61 nm, respectively. Thus, if the 
separation between an atom and a plate is a = 3 nm, then the plate of d = 8 nm thickness 
can be already considered with a good accuracy as a semispace. 



VI. CALCULATION OF THE VAN DER WAALS INTERACTION FOR 
HYDROGEN ATOM OR MOLECULE EXTERNAL TO MULTIWALL CARBON 
NANOTUBE 



The multiwall carbon nanotube can be modelled by a graphite cylindrical shell of some 
length L, external radius R <^ L and thickness d < R. In doing so the hexagonal layers 
of graphite crystal lattice form the external surface of a cylinder and the internal sections 
concentric to it. The crystal optical axis z is perpendicular to the surface of the cylinder 
at each point. The above derived Lifshitz-type formula (fT^ is applicable to the case of 
multiwall carbon nanotube if its thickness d is large enough (typically D > 3nm), so that 
the nanotube contains sufficiently many layers. Then it is possible to neglect the atomic 
structure of graphite and to describe it in terms of dielectric permittivity. 

For convenience in numerical computations we rewrite Eq. (fT^ in terms of dimensionless 
variables ()30p representing the free energy of the van der Waals interaction with a cylinder 
in the form 

F^Eia,T) = (35) 

where 



C3^(a,T) 



R 



AR + 3a 



8 \ R + a \2{R + a) 



a(0) 



poo 

+ ^ct{iCi^^c) / dyye' 
1=1 -^ii 



2{R + a) 



X 



(36) 



The reflection coefficients were defined in Eq. ()32j) (with index s related to the case of a solid 
cylinder) and in Eq. ()34|) (with index p related to the case of a cylindrical shell of thickness 
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d). 

Let us first compare the van der Waals interaction between hydrogen atom or molecule 
with a graphite semispace and a solid cylinder. The differences of the interaction strength 
with a semispace and a cylinder can be characterized by a parameter 6 = (C| — C|)/C|. 
A few results for a graphite cylinder with R = 50 nm, calculated by Eqs. (jSHI), PTjl- (PT|). 
()22|) and dielectric permittivities given by the solid lines of Fig. 2, are presented in Table II 
(columns 2-4 and 5-7 are related to the cases of hydrogen atom and molecule, respectively). 
As is seen from Table II, at short separations of about a few nanometers there are only minor 
differences between C| and C3. With increase of a, however, the magnitude of 6 quickly 
increases. This takes place for both hydrogen atom and molecule. 

It is interesting to follow the dependence of the van der Waals coefficient C| on R for atoms 
and molecules located at different separations from the cylinder surface. These computations 
were performed with Eqs. (jSHl), ®! (EH); and the same data for graphite dielectric 
permittivities. The results are presented in Fig. 5a (for hydrogen atom) and Fig. 5b (for 
hydrogen molecule) where the lines 1, 2 and 3 are pictured for separations a = 3nm, 5nm 
and 10 nm, respectively. It is seen that with the increase of R the van der Waals coefficients 
are also increasing. 

Now consider the cylindrical shell of radius R and thickness d with the longitudinal cavity 
of a radius R — d. This is evidently a better model for a multiwall carbon nanotube. In 
Fig. 6 we present the computation results for the interaction between a hydrogen atom 
and a cylindrical envelope with R = 20 nm as a function of envelope thickness d (atom is 
located at a separation a = 5nm from the external surface of the cylindrical shell). The 
computations were performed by Eq. (jHBj) using the same procedure as above. The value 
d = 20 nm corresponds to the case of a solid cylinder. It is interesting, however, that 
already at c? = 11 nm the magnitude of C| is only 1% lower than the one obtained for the 
solid cylinder of R = 20 nm radius. For less thickness of the cylindrical shell the smaller 
values of the van der Waals coefficient are obtained (the same is true also for a hydrogen 
molecule). Note that we do not extend the line of Fig. 6 for thicknesses less than 3nm 
where the macroscopic description of graphite in terms of dielectric permittivity may be not 
applicable. 
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VII. COMPARISON BETWEEN THE FREE ENERGIES OF HYDROGEN 
ATOMS INSIDE AND OUTSIDE OF MULTIWALL CARBON NANOTUBES 



The obtained above Lifshitz-type formulas (fT^. provides a good approximate de- 
scription of the van der Waals interaction when a microparticle is located outside of a 
cylindrical shell. Let us now consider a microparticle inside of the same shell. In this case 
the van der Waals free energy can be approximately calculated by the method of pairwise 
summation of the interatomic potentials with subsequent normalization of the obtained in- 



teraction coefficient using the known case of microparticle near a semispace 







4^. For a 



microparticle outside of an arbitrary macrobody v this method leads to the expression 

Fl(a,T)^-'-!^^l%, (37) 

where r is the separation between the microparticle and an atom (molecule) of the macro- 
body. 

To determine the accuracy of Eq. (|H7j). let us apply it in the case of hydrogen atom 
outside of a solid graphite cylinder at a separation a [to which Eq. (jHUj) is also applicable]. 
Then Eq. (IH7j) is rewritten as 

F, =F,(a.T)« —J^ dej^ dzj^^^^^ (38) 

where sin^^ = R/ {R + a), R is the cylinder radius, and pi,2(^) are the two solutions of the 
equation 

+ (i? + af - 2p{R + a) cos 9 = R^. (39) 
After the integration over z and p Eq. ()H8|) takes the form 

'1 1 ■ 

The numerical computations by Eq. ()40j) demonstrate that for a cylinder with R = 50 nm 
the results, obtained by the method of additive summation, differ by less than 1% from the 
results, obtained by the Lifshitz-type Eq. (|35p. within the separation range a < 8nm. At 
a = 10 nm the free energies computed by the two formulas differ for 1.35%, and at a = 50 nm 
by 16%. Hence the method of additive summation works well at small separations between 
an atom and a cylindrical surface. This makes it reasonable to apply this method for 
hydrogen atom inside of a multiwall carbon nanotube. 
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o r9m 

Fr\a,T) ^ --C',{a,T) dO 



(40) 



We consider a hydrogen atom inside of a nanotube with thickness d and internal radius 
Rq = R — d aX a, separation a from the internal surface. In accordance with Eq. (pTf). the 
free energy of the van der Waals interaction is 



pria, T) . I ,e I J^^^^^ (41) 



where the integration limits are given by 



pi(0) = -(i?o_a) cose+ jRl-{Ro-aysmH 



p^[e) = ~{Rq- a) cose+^ {Ro + df - (i?o - a)2 sin^ 9. (42) 
After the integration over z and p Eq. ()4H1 leads to 



de 



1 1 



(43) 



In Fig. 7 we present the results of numerical computations by Eq. ()43|) for the hydrogen 
atom inside of the hypothetical nanotube with the internal radius Rq = 10 nm and external 
radius R = 50 nm. The free energy of the atom-nanotube interaction is plotted in Fig. 7 as a 
function of atom position between the opposite points of the internal cylindrical surface. The 
atom positions closer than 3 nm to the internal surface are not reflected in the figure (their 
consideration would demand a more exact treatment of the atomic structure of graphite). As 
is seen from Fig. 7, the free energy reaches a maximum on the cylinder axis, where the van der 
Waals force acting on an atom is equal to zero in accordance with symmetry considerations. 
This equilibrium state is, however, unstable and under the influence of fluctuations the 
hydrogen atom will move to positions with lower free energy near the internal cylindrical 
surface of a nanotube. 

Now we are in a position to compare the free energies of hydrogen atoms located outside 
and inside a multiwall carbon nanotube in order to decide which position is preferable 
energetically. In Fig. 8 the calculation results for the differences of free energies F^^^^ and 
F^'*"'* are presented as a function of thickness of the nanotube. In doing so we consider 
both atoms, internal and external, situated at a separation a = 3 nm from the internal and 
external surfaces of a nanotube, respectively. The solid line in Fig. 8 is related to the fixed 
internal radius of the nanotube Rq = 10 nm, and in this case the external radius increases 
together with thickness of the nanotube d. The dashed line is for a fixed external radius 
R = 50 nm and decreasing internal radius with the increase of d. The computations were 
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performed with Eq. ()43p for a position of the atom inside the nanotube and with Eq. ()35|) 
for position of the atom outside the nanotube. 

As is seen from Fig. 8, in all cases the difference between the external and internal free 
energies of the van der Waals interaction is positive. What this means is the position of a 
hydrogen atom inside a multiwall carbon nanotube is preferable energetically. Comparing 
the solid and dashed lines in Fig. 8, we conclude that for nanotubes of fixed thickness d the 
potential well for the hydrogen atom inside a nanotube is deeper if nanotube has a smaller 
external radius R. This is an encouraging result which points to the possibility of hydrogen 
storage inside carbon nanostructures. 

VIII. CONCLUSIONS AND DISCUSSION 

In the above we have widened the scope of the Lifshitz theory of the van der Waals force 
by considering new configurations of much interest which have not been explored previously. 
The first to be investigated was the van der Waals force between an atom or molecule 
and a plane surface of an uniaxial crystal perpendicular to the crystal optic axis. For this 
configuration the exact expression for the free energy of the van der Waals and Casimir- 
Polder interaction is given by Eq. (P) with the reflection coefficients (0) (for the case of a 
microparticle near a semispace) or (for a microparticle near a plate of finite thickness). 
We next derive the approximate Lifshitz-type formula for the free energy of the van 
der Waals interaction between microparticle and solid cylinder or cylindrical shell having 
a longitudinal concentric cavity. This cylinder may be made of isotropic material or of an 
uniaxial crystal. The accuracy of the obtained formula was shown to be of about 1% at 
microparticle-cylinder separations less than one half of a cylinder radius. 

The above extensions of the Lifshitz formula for microparticle-wall interaction were ap- 
plied to the case of hydrogen atom or molecule near a graphite surface. For this purpose 
the dielectric permittivities of graphite along the imaginary frequency axis were found by 
the use of tabulated optical data for the complex refractive index. In doing so different 
sets of data were analyzed and necessary extrapolations to high and low frequencies were 
done. Together with the use of hydrogen atomic and molecular dynamic polarizabilities, this 
allowed us to calculate the van der Waals interaction between hydrogen atom or molecule 
and graphite semispace, graphite flat plate of finite thickness or solid graphite cylinder and 
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cylindrical shell. In particular, the influence of the thickness of the plate on the van der 
Waals interaction was investigated. 

The calculation results for the atom-cylinder case were used to model the van der Waals 
interaction between hydrogen atoms or molecules and multiwall carbon nanotube with suffi- 
ciently large number of layers. In particular, the dependence of the van der Waals interaction 
of the atom-nanotube case on nanotube thickness was investigated. Notice that the devel- 
oped formalism is not applicable to single- or twowall nanotubes where the atomic structure 
of the wall should be taken into account. In this case the van der Waals force can be 
computed in the framework of density functional theory {45!. I46. 42 1. 

Finally, we have compared the free energies of the van der Waals interaction between a 
hydrogen atom and multiwall carbon nanotube for the cases when atom is located outside 
or inside of the nanotube. It was shown that atoms situated inside of a multiwall nanotube 
possess lower free energy in a wide region of nanotube thicknesses, i.e., such a position is 
energetically preferable. This conclusion is promising for the possibility of using carbon 
nanotubes for the purpose of hydrogen storage. 

Many other opportunities for application of the obtained generalizations of the Lifshitz 
formula in physics of dispersion forces are possible. 
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FIG. 1: Schematic of the cylinder of radius R made of a uniaxial crystal and having a longitudinal 
concentric cavity of radius R — d. This cylinder is concentrically placed into a cylindrical cavity of 
radius R + a m. the infinite space filled with an isotropic substance. 
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FIG. 2: Dielectric permittivity of graphite along the imaginary frequency axis in (a) the hexagonal 
layer and (b) perpendicular to it, as a function of frequency. Solid and dashed lines in (b) are 
obtained with the optical data of Ref. [3] and Ref. respectively. 
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FIG. 3: Dependence of the van der Waals coefficient C| on separation of (a) hydrogen atom and 
(b) molecule, from graphite semispace. The solid and dashed lines are obtained with the optical 



data of Ref. 



and Ref. 



4l| |. respectively. 
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FIG. 4: The ratios of the free energies for the van der Waals atom-plate to atom-semispace 
interaction as a function of plate thickness for hydrogen atom located at different separations from 
the graphite surface (lines 1, 2, 3 and 4 are for separations a = 3nm, 10 nm, 20 nm and 50 nm, 
respectively) . 
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FIG. 5: Dependence of the van dcr Waals coefficient C| on the cyhnder radius for (a) hydrogen 
atom and (b) molecule, located at different separations from the graphite cylinder (lines 1, 2 and 
3 are for separations a = 3nm, 5nm and 10 nm, respectively). 
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FIG. 6: Dependence of the van der Waals coefficient C| on thickness of the cyhndrical sheU with 
an external radius i? = 20 nm for hydrogen atom at a separation a = 5 nm from the shell. 
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FIG. 7: The van der Waals free energy for hydrogen atom inside of the carbon nanotube with 
internal radius = 10 nm and external radius R = 50 nm as a function of the atom position 
between the opposite points of the internal cylindrical surface. 
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FIG. 8: Difference of the free energies of hydrogen atoms situated outside and inside of the 
multiwall carbon nanotube as a function of nanotube thickness. The sohd and dashed hnes are 
for the nanotubes with a fixed internal radius Rq = 10 nm and fixed external radius R = 50 nm, 
respectively. 
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TABLE I: The values of strengths and eigenenergies of oscillators for hydrogen atom in the frame- 
work of lh(; 10-oscilbil(ir iiKKiel. 



j 


9j 


UJaj (a.c.) 


1 


0.41619993 


0.37500006 


2 


0.08803654 


0.44533064 


3 


0.08993244 


0.48877611 


4 


0.10723836 


0.56134416 


5 


0.10489786 


0.68364018 


6 


0.08700329 


0.89169023 


7 


0.06013601 


1.2698693 


8 


0.03259492 


2.0478339 


9 


0.01199044 


4.0423429 


10 


0.00197021 


12.194172 
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TABLE II: Magnitudes of the van der Waals coefficients and C3 and their relative differences S 
(see text) for the interaction of hydrogen atom or molecule with a graphite semispace or a cylinder 
with radius i? = 50 nm. 



a 


H 


H2 


(nm) 


C|(a.u.) 


C|(a.u.) 


6{%) 


C|(a.u.) 


C|(a.u.) 


6{%) 


3 


0.09882 


0.09471 


4.2 


0.1317 


0.1262 


4.2 


5 


0.09416 


0.08792 


6.6 


0.1248 


0.1166 


6.6 


10 


0.08316 


0.07322 


12.0 


0.1088 


0.09584 


11.9 


20 


0.06652 


0.05301 


20.3 


0.08526 


0.06801 


20.2 


30 


0.05516 


0.04047 


26.6 


0.06970 


0.05118 


26.6 


40 


0.04704 


0.03214 


31.7 


0.05885 


0.04025 


31.6 


50 


0.04098 


0.02631 


35.8 


0.05090 


0.03270 


35.8 
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